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Cn . We study flavor-conserving radiative decays of double heavy baryons using a manifestly Lorentz 

covariant constituent three-quark model. Decay rates are calculated and compared to each other in 
the full theory, keeping masses finite, and also in the heavy quark limit. We discuss in some detail 
hyperfine mixing effects. 
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I. INTRODUCTION 

A first observation of the double charmed baryon S+ (3519) by the SELEX Collaboration at Fermilab 1] stimulated 
theoretical studies of double heavy baryons (DHBs). Up to now the study of DHBs has mainly focussed on their mass 
spectra and their semileptonic decays (for an overview see e.g. Refs. [3,|3|). In particular, significant progress has been 
achieved in the analysis of the DHB semileptonic weak decays. The current-induced flavor-changing double-heavy 
baryon transitions have been analyzed in a number of model app roaches. These include effective field theories based 
on heavy quark spin symmetry [Ji7|], three-quark models |8l-lll| , quark-diquark models |12l . |l3| , and nonrelativistic 
QCD sum rules 0, llJ]- Recently Ql we have presented a comprehensive analysis of the semileptonic decays of 
DHBs using a manifestly Lorentz covariant field theory approach termed the relativistic constituent three-quark 
model (RTQM) |^, 8., 2M ■ We considered all possible current-induced spin transitions between double-heavy baryons 
containing both types of light quarks - nonstrange q — u,d and strange s. These involved the flavor-changing 
transitions be — > cc and 66 — ^ be. Form factors and decay rates have been calculated and have been compared 
to each other in the full theory with all masses finite, and also in the heavy quark limit (HQL). Such an analysis 
is important because the semileptonic decays of DHBs provide yet another opportunity to measure the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element Vet- This is particularly true since the transition matrix elements between 
double-heavy baryons obey spin symmetry relations in the heavy quark limit in addition to a model independent zero 
recoil normalization of the relevant transition matrix elements. 

In this paper we continue the study of DHB properties in the RTQM 0, [I^. In particular, we analyze flavor- 
conserving radiative transitions between ground state DHBs: 1/2+ — > 1/2+ and 3/2+ — > 1/2+. The first estimate of 
the radiative decay widths of the DHBs in the heavy quark limit including hyperfine mixing effects has been done in 
Ref. ^] (for radiative transitions between DHBs see also the early paper [IB])- A detailed analysis of DHB decays 
containing only the (6c) heavy quark configuration has been considered recently in Ref. [I3|. As in our recent paper [J], 
we take the DHBs to be bound states of a light quark and a double-heavy (Q1Q2) diquark. 

The origin of the hyperfine mixing for double heavy baryons is the one-gluon exchange interaction between the 
light and heavy quarks in the DHB states containing two different heavy quarks — b and c. It leads to mixing of the 
states containing spin-0 and spin-1 heavy quark configurations. As shown in Refs. |10l . I 111 . Il7l | hyperfine mixing has 
a big impact on the decay properties of double heavy baryons. Both the weak semileptonic and the electromagnetic 
decay rates involving mixed DHB states are significantly enhanced or reduced relative to the rates involving unmixed 
states. 

The RTQM can be viewed as an effective quantum field theory approach based on an interaction Lagrangian of 
hadrons interacting with their constituent quarks. From such an approach one can derive universal and reliable predic- 
tions for exclusive processes involving both mesons composed of a quark and antiquark and baryons composed of three 
quarks. The coupling strength of a hadron H to its constituent quarks is determined by the compositeness condition 
Zh = [la. Il9l| where Zh is the wave function renormalization constant of the hadron H. The quantity Z^j is the 
matrix element between the physical particle state and the corresponding bare state. The compositeness condition 
Zh = enables one to represent a bound state by introducing a hadronic field interacting with its constituents so that 
the renormalization factor is equal to zero. This does not mean that we can solve the QCD bound state equations but 
we are able to show that the condition Zh — provides an effective and self-consistent way to describe the coupling 
of a hadron to its constituents. One starts with an effective interaction Lagrangian written down in terms of quark 
and hadron variables. Then, by using Feynman rules, the S'-matrix elements describing hadron-hadron interactions 
are given in terms of a set of quark level Feynman diagrams. In particular, the compositeness condition enables one 
to avoid the problem of double counting of quark and hadronic degrees of freedom. The approach is self-consistent 
and all calculations of physical observables are straight-forward. There is a small set of model parameters: the values 
of the constituent quark masses and the scale parameters that define the size of the distribution of the constituent 
quarks inside a given hadron. 

The main objective of the present paper is to present an analysis of all possible electromagnetic transitions between 
ground state DHBs containing both types of light quarks - nonstrange q ^ u,d and strange s. The paper is structured 
as follows. First, in Sec. II we review our relativistic constituent three-quark model approach (for more details see 
6-g- |3|) including a discussion on how to obtain a gauge invariant coupling of the photon in our model. In Sec. Ill 
we discuss in more detail various aspects of the radiative decays of DHBs. We discuss the calculation of the relevant 
radiative transition matrix elements and analyze the consequences of taking the heavy quark limit for the radiative 
transitions. In Sec. IV we discuss in some detail hyperfine mixing effects in the radiative decays of DHBs. Sec.V 
contains our numerical results which are compared to the predictions of a naive nonrelativistic quark model that has 
the same spin-flavor symmetry group as our DHB currents in the nonrelativistic limit. We also compare the results 
of the full finite mass calculation with results derived in the HQL. In addition we compare our results for 6c — > be 
radiative transitions with recent quark model results [iTJ. Finally, in Sec. VI we present a brief summary of our 
results. 



II. FRAMEWORK 
A. Lagrangian 

For the evaluation of the radiative decays of DHBs we will consistently employ the RTQM [3, [l^ . The model is 
based on an interaction Lagrangian describing the coupling between a baryon -B(gi(j2<Z3) and its constituent quarks 
Qi, q2 and 93- For J^ = ^ and | baryons the Lagrangians read: 



'CLt(3^) = 9bB{x) dxi dx2 dx3F{x,Xi,X2,X3)JBixi,X2,X3) 

+ 9b-B*^{x) dxidx2dx3F{x,xi,X2,X3)J'^,{xi,X2,X3)+h.c. (1) 

where Jb and Jb* are interpolating three-quark currents with the quantum numbers of the relevant baryon B{^ ) 

and i?*(| ). Note that the spin 3/2 spinor corresponding to the i3*(| ) field satisfies subsidiary Rarita-Schwinger 
conditions (see further details in Appendix A). 
One has 

Jb{xi,x2,x3) = e''''''"'Tiq''^{x3)Qt'{xi)Cr2QTix2), (2) 

J^,{x,,x2,X3) = e''''''''-^r,q'''{x3)Q'','ix,)CTt^Q^^'{x2), (3) 

where the ri_2 are strings of Dirac matrices, C is the charge conjugation matrix C = j'^j'^ and the at (i=l,2,3) are 
color indices. F{x, a;i, 0:2, 2^3) is a nonlocal scalar vertex function which characterizes the finite size of the baryons. 
The full Lagrangian 

C,^,) ^ /:f_(x) + cZ^'\:c) + £:r°'"^''(^) (4) 

needed for the calculation of the radiative decays of DHBs includes the free parts of the baryons and the constituent 
quarks 

U,,{x) ^ B{x)VbB{x) ~ Bl{x)V%^,Bl + ^ i^{x)V^,ij{x), (5) 

where 

V^, = i ^ - m^ , 
Vb = i ^ - TUB , 
V^", = g^"'{i ^~mB-)~ ih^d" + -/"df") + Y^ ^Y " ^5.7^7" . (6) 

The baryon and constituent quark masses are denoted by mB{B*) and m^, respectively. 
The electromagnetic interaction Lagrangian contains two pieces given by 

/•em _ ^cm(l) ^cm(2) /-.x 

which are generated after the inclusion of photons. The first term >C°™ is generated via minimal substitution in the 
free Lagrangian £frco: 

9^^ ^ [di" - ie^A")^ , d"'^ -^ {d" + ie^A")'^ , (8) 

where ^ stands for B, B* , q, Q, and e* is the electric charge of the field ^. The interaction Lagrangian /3™ reads 



^cni(l) 



{x)^eBB{x)^B{x)~eB'Bl{x)Ug^'''^r4Y-l''A'' -YA''\bI{x)+ ^ e^^:{x)^i:{x) . (9) 

The second electromagnetic interaction term C^^ is generated when one gauges the nonlocal Lagrangian Eq. ([T|). 
The gauging proceeds in a way suggested and extensively used in Refs. [H, [20, [2l|. In order to guarantee local 



electromagnetic gauge invariance of the strong interaction Lagrangian one multiplies each quark field in £f*^ with a 
gauge field exponential. One then has 

X e-^'^Qi^^^i'^'-P) Qi'(xi)Cr2 e-*'^«2-f(^2,:.,P) g^s^^.^) ^ 

+ gB'Blix) /"dxi/"da;2/"dx3F(x,xi,X2,X3)e'^i"^"-^rie-*^'^("^'"'^)(7''^(x3) 

where 

Iix^,X,P)= f dz^Af'iz). (11) 

X 

An expansion of the gauge exponential up to a certain power of A^^ leads to the terms contained in C°^ . 

The full Lagrangian consistently generates all the required matrix elements of the radiative decays of the DHBs. 
The relevant transitions can be represented by a set of quark loop diagrams. In the evaluation of the quark loops we 
use the free fermion propagator for the constituent quarks as dictated by the free quark Lagrangian discussed above. 
One has 

zS^{x-y) = {0\T^ix) ^iy)\0) = J ^ e-^'f^-^) S^{k) (12) 

where 

SAk) = V— (13) 

is the usual free fermion propagator in momentum space. We avoid the appearance of unphysical imaginary parts 
in Feynman diagrams by postulating that the baryon mass is less than the sum of the constituent quark masses 
"^B(s*) < '^g + "^Qi + ^f^Q2 which is satisfied in our calculation. We mention that we have recently introduced a 
further refinement of our model in that we can now include quark confinement effects [22] ■ 
The free propagators of the baryon fields in momentum space are given by 

SB{k) = \^, (14) 



S%' (k) - ^— -9"" + o7^7^ + ,r^ + -^, ■ (15) 



Next we consider in detail the required building blocks of the strong interaction Lagrangian >C?*j — the vertex function 
F, the interpolating three-quark currents Jg and Jg, , and the baryon-quark coupling constants g^ and g^* ■ 

B. Vertex function 

The vertex function Fg is related to the scalar part of the Bethe-Salpeter amplitude and characterizes the finite 
size of the baryon. In our approach we use a specific form for the vertex function given by 






F(a;,a;i,a;2,a;3) = N 6^^> {x ~ } ^w,x,) '^l }^{x^ - x^Y ] , (16) 

i<j 

which is Poincare-invariant. $ is a nonlocal correlation function involving the three constituent quarks with masses 
mi, 777.2, TT13', A'^ = 9 is a normalization factor. The variable 77;^ is defined by Wi = 777^/(7771 + m,2 + 7^73). 



The Fourier transform of the correlation function $ I J2{^i~ ^jY I can be calculated by using Jacobi coordinates. 

\i<0 J 

One has 

= {^-.fs^'^v-Y^v^H-il^iD. (17) 

where the Jacobi coordinates are defined by 

^1 " ^ + 75 ^1^3 - ^ 6(2w2 + w-i) , 
^2 " ^ + 75 ^1^3 + ^ 6(2wi + W3) , 
2:3 = a; - ^ ^i{wi + W2) + ^ 6(^1 - u'2) • (18) 

The corresponding Jacobi momenta read 

P = P1+P2+P3, 

h = ■^■W3{pi+P2)~ ^{Wi+W2)p3, 

h = -^i'2^i'2+W3)pi + ^{2wi+W3)p2 + ^iwi-W2)p3, (19) 

3 / 

^2 



where, according to Eq. p^ . ^ WiXi — x. Since the function $( J2i^i ~ ^j)^ ) i^ invariant under translations its 

1=1 \i<j / 

Fourier transform only depends on two four-momenta. The function $(— ^^ — ^|) in Eq. (J17p will be modelled by 
a Gaussian form in our approach. The minus sign in the argument is chosen to emphasize that we are working in 
Minkowski space. Our choice is the Gaussian form 

$(-/? _ q) ^ exp (l8{ll + q)/A') (20) 



where the parameter A characterizes the size of the DHB. Since If and I2 turn into —if and — /^ in Euclidean space 
the form (j20|) has the appropriate fall-off behavior in the Euclidean region. 

C. Three— quark currents 

In the so-called Qi(52~basis the DHBs are classified by the set of quantum numbers {J^,Sd), where J^ is the 
spin-parity of the baryon state and Sd is the spin of the heavy diquark. There are two types of heavy diquarks - 
those with Sd — (antisymmetric spin configuration [Q1Q2]) and those with 5"^ = 1 (symmetric spin configuration 
{QiQ2})- Accordingly there are two J^ = 1/2+ DHB states. We follow the standard convention and attach a prime 
to the Sd — states whereas the Sd = i states are unprimed. The J^ — 3/2 "*" states are in the symmetric heavy 
quark spin configuration. In Table I we list the quantum numbers of the double-heavy baryons including their mass 
spectrum as calculated in [l3j . 

As we have discussed in our recent paper Q, there is a mass inversion in the (1/2+) mixed flavor states (S^c, 2[,^) and 
(fifjc, ^bc) in that Af (SJ^^) > M(Sftc) and M{Wi^^) > M(rif,c) even though the heavy diquarks satisfy the conventional 
hyperfine splitting pattern m((,c)s=i > 'n^(bc)s=o- This inversion is a feature of all models that have attempted to 
calculate the mass spectrum of double-heavy baryons @, M, SS l24l - [2a |. In particular, the inverted mass hierarchy 
implies that one can only expect substantial flavor-changing branching ratios for the two lowest lying states S^c 
and flbc whereas the rates of the higher lying states SJ,^, S^^ and fi^,^, fll^ will be dominated by flavor-preserving 
one-photon transitions to the lowest-lying states She and i^bc- One of the purposes of the present paper is to analyze 
the strength of the one-photon transitions between the Sd = and Sd = I double-heavy baryon states. In the HQL, 
the photon couples to the light quark only, and therefore one-photon transitions between the Sd = and 5^ = 1 
double-heavy baryon states such as S^^ -^- Ebc + 7 are forbidden in this limit. For finite heavy quark masses one- 
photon transitions between the Sd = and Sd = ^ double-heavy baryon states will occur at a somewhat reduced rate 
which, however, is very likely to exceed the flavor-changing weak decay rates of these states |16j . 



Following the suggestion of Ref. [10] (see also discussion in [Tl|) we also consider hyperfine 2^^ — 5bc and J7J,^ — flbc 
mixing induced by one-gluon exchange interactions. We define the mixed states through the unmixed states using a 
unitary transformation [l3, llll with the mixing angles 9s or 9n : 

(Bt\f cosOb sin 0s \ / i?L \ 

\Hc/ \ - sin^B cos^B J \BbcJ ' 

where B = S or n. We treat the mixing angle 9b as a quantity of order 0{as) where as is the QCD coupling constant. 
By B'^ and S' we denote the mixed states with the mass hierarchy to^, < ttlb < itib' < i^Bh ■ The masses of the 
mixed states to^^ and vriBi differ from the masses of the unmixed states by small hyperfine splitting corrections. They 
are expressed through the unmixed masses and mixing angle as: 



TOBh = fnB' + {rriB' - rnB) 



sin^9 
cos 261 ' 
sm^9 



rUBi = mB-{mB'-mB) xt; ■ (22) 

cos 29 

When diagonalizing the mass matrix of the unmixed states the diagonal elements are driven apart such that the mass 
difference of the mixed states is larger than that of the unmixed states. One has 

"^^^ "'^' = cos 20 ■ ^2^) 

This leads to a further enhancement of the widths of the radiative transitions Bj^^ — > Bj^^ + 7 between mixed states, 
because the photon transition rate is proportional to (Mi — M2)^ (see Appendix A). 

The mixing angle 9b corresponds to the combination 9 + 30° in [ll|: 9b = 9 + 30°. Here 9 is the angle that 
rotates the mixed B^, B^ states into the B, B' states — bound states of 5-quark and heavy-light cq diquark (so-called 
qc-basis). The angle 30° corresponds to a further rotation of the B,B' states into the unmixed B'^B states defined 
in the 6c-basis (i.e. bound states of light quark and heavy be diquark). In the quark model calculation of |ll| one 
obtains sin^s = 0.431 {9-e. = 25.5°) and sin0Q = 0.437 (0a = 25.9°). Using these values of the mixing angle 9 and the 
masses of unmixed states we deduce the following values for the masses of the mixed states: 

m^H = 6.972 MeV , m^i = 6.924 MeV , 

"ba "be 

m^b ==7.125 MeV, m^i^ = 7.079 MeV . (24) 

The interpolating currents of the DHB states BqQ^q^ are constructed in the form of a light quark q°'^ coupled to a 
heavy diquark D'q q . One obtains 

J,Q.Q. = r^.Q.93^^QWv D'^q\q.. - ^'^^'^^'^^ {QrCT^.Q..QT) ■ (25) 

For the (5 ,0), (i ,1) and (| ,1) states we use the simplest currents — the pseudoscalar J^ , the vector J^ and 



(26a) 
(26b) 
(26c) 

(27a) 
(27b) 

(27c) 

where '>Pq,Qi,Q2 ^^^ ^^e upper components of the Dirac quark spinors and the ai are the Pauli spin matrices. Note 
that the spin-flavor wave function coincides with the nonrelativistic limit in the HQL. In the nonrelativistic limit our 



^ t^LliJ-CliUO, i (^OJ-*t;\jlJl V V^iJ . 




jP _ 010203 

'^qQlQ2 - ^ 


q'''{QrCj,QT), 


7V .010203 

-'gQiQ2 - ^ 


ri'q'''{QVCj^Q' 


J^V 010203 


q^'{Qrc-f^QT)- 


n the heavy quark limit the above currents reduce to 




aP _ 010203 

J9Q1Q2 - ^ 


rqUrQ\^2rQ\), 


aV _ 010203 

hQiQi - ^ 


ar^ (V^^asaVgV 


1-V _ 010203 


<3(^^.^,^^-), 



DHB currents have a one-to-one correspondence to the naive quark model baryon spin-flavor functions (up to overall 
factors) which are displayed in Table II. Further details on the naive quark model and how to evaluate the radiative 
transition amplitudes in this framework can be found in Appendix B. 

According to the definition (PT|) the interpolating currents of the mixed B^ and S' states are given by 



cos 9 b sin 6 b 





(28) 
-sin 6*5 cos 6*5 / V, gBJqbc 

D. Baryon-quark coupling constant 

The coupling constants g^ (■0 = B,B*) are determined by the compositeness condition [2, UM, llM UM- The 
compositeness condition implies that the renormalization constant of the hadron wave function is set equal to zero, 
i.e. 

Z^, = I - gli:'^im^) = . (29) 

S^(m^) is the derivative of the baryon mass operator shown in Fig.l. 

In case of the 3/2+ states the function E^. (p) is subtracted from the nonvanishing part of the mass operator S^'C {p) 
proportional to the Minkowski metric tensor g'^": 

Ui_,{p,s*)T.'^g,(p)uy{p,s*) = Uf_,(p,s*)'EB'ip)g'^''u^ip,s*) , (30) 

where u^{p, s*) is the |-spinor. Other possible Lorentz structure on the r.h.s. of Eq. (j30[) vanish due to the Rarita- 
Schwinger conditions. Note that the compositeness condition is equivalent to a Ward identity relating the electro- 
magnetic vertex function at zero momentum transfer to the derivative of the mass operator (see details e.g. in [3|). 
Explicit expressions for the baryon mass operators are given in Appendix C. 

III. RADIATIVE DECAYS OF DOUBLE HEAVY BARYONS 

A. Matrix elements 

In our approach the radiative decays of DHBs are described by the set of Feynman diagrams shown in Fig. 2. The 
three "triangle" diagrams Figs.2(a)-2(c) are generated by the coupling of the constituent quarks with the photon. The 
two "bubble" diagrams in Figs. 2(d) and 2(e) are generated by gauging the nonlocal strong Lagrangian (see discussion 
in Sec. II). Finally the two "pole" diagrams in Figs. 2(f) and 2(g) are generated by the direct coupling of the initial/final 
baryon with the photon. Due to our explicit construction the DHB radiative matrix elements are explicitly gauge- 
invariant. The pole diagrams vanish for the radiative transitions (^ ,0) — )■ (^ , 1) and (| ,!)—)• {\ ,0) due to the 
orthogonality of the heavy diquark spin wave functions. For the same reason the photon does not couple to the light 
quarks in these modes implying that the corresponding triangle and bubble graphs vanish. For the (| , 1) — > (^ ,1) 
transitions the pole diagram in Fig. 2(g) vanishes due to the Rarita-Schwinger conditions for the |-spinor. All these 
statements are true for photon transitions between unmixed states. Photon transitions between mixed DHB states 
will be discussed in section IV. 

We continue with a summary of some useful analytical results. The on-shell matrix elements for the radiative 

UB2 {P2,S2) Af_,{pi,P2) UBi ipi,si) , 
UB2{P2,S2)k^^{pi,P2)uB.{pi,S*), (31) 

where ub{p, s) and u^* {p, s*) are the spin ^ and | spinors with the normalization (see further details in Appendix A): 

ub{p-,s)ub{p,s) = 2mB , u'^,{p,s*)uB*ip,s*) = -2mB' ■ (32) 

The momenta of the final state photon, the initial and final state baryon are denoted by q, pi and p2, respectively, 
where q — pi ^ P2 and where s and s* are spin indices. Due to gauge invariance the electromagnetic vertex function 
A^(pi,P2) is orthogonal to the photon momentum q^K^(pi^p2) = 0. As a result the vertex function A^(pi,p2) is 
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8 

given by the sum of the gauge- invariant pieces of the triangle (A), the bubble (bub) and the pole (pol) diagrams, 
while the nongauge-invariant parts of these diagrams cancel in the sum: 

A^{puP2) = A^^^(pi,p2) + A^,3^b(pi,P2) + A^p„i(pi,P2) (33) 

The contribution of each diagram can be split into gauge invariant and gauge variant pieces by introducing the 
decomposition 

such that 7^ q^ — 0, pj^ q^ — 0, and gj^^^ q^ — 0, where pi is pi or p2- The vertex function A^(pi,p2) can then be 
expressed in terms of j^, the p,^ and g^^. Note that all matrix elements are finite for real photons {q^ = 0). Doing 
our calculations we start with q^ ^ and then take the limit (7^ — > 0. Explicit expressions of the electromagnetic 
vertex functions can be found in Appendix C. 

B. Heavy Quark Limit 

In the HQL the masses of the heavy quarks are taken to infinity (ttiq — >■ 00). In this limit the spins of the double- 
heavy diquark and the light quark in the DHB states decouple leading to a much simplified transition structure. In 
particular, the transitions amplitudes (^ ,0) — > (^ ,1) and (| , 1) — > (5 ,0) vanish as ©(l/mg) since the photon 
coupling to the heavy quarks involves a spin-flip factor proportional to the magnetic moment of the heavy quark 
given by fiq = cq / (2mQ) . This is in full agreement with the nonrelativistic quark model (for more details see the 
discussion in Appendix B) where the corresponding amplitudes are found to be proportional to the difference of the 
magnetic moments of the heavy quarks /ic — /ift. The transition amplitude (| , 1) — > (^ ,1) survives in the HQL since 
the photon can now couple to the light quark. Again, this is consistent with the nonrelativistic quark model. The 
structure of the amplitude of the (| , 1) — ?• (^ ,1) transition significantly simplifies in the HQL. Only the triangle 
diagram in Fig. 2(a) contributes to the transition amplitude since Fig. 2(a) represents the direct coupling of the light 
quark with the photon. The contribution of Fig. 2(a) scales as 0(1) in the inverse heavy quark mass expansion. The 
other diagrams are suppressed in the HQL. In particular, the triangle diagrams in Figs. 2(b) and 2(c) contribute only 
at ©(I/toq) since they represent direct couplings of the heavy quarks with the photon. The same holds true for the 
bubble diagrams Figs. 2(d) and 2(e) involving a nonlocal photon-light quark coupling. 

Following ideas developed in our paper ^3], we choose the momenta of the initial and final DHB as p^ = {niQ-^ + 
mQ^)v^ and Pj = ('^Qi + '"^Qi)'^'^ — ("^Qi + ^^Q2)^^ + '"^ where r is a small residual momentum in the sense that 
r^ ~ 0{1) when mq^ — >■ cxd. with these assumptions the heavy quark propagators simplify in the HQL. One has 

SQ^(h±Pnr) -> 't'^^j-Ze ' ^^^) 

where -qi = rnqJimQ^ + toqJ. 

In the HQL and at q^ — the explicitly gauge invariant transition amplitude (| ,1)— 7>(i ,l)is given by 

A^. - r^i /a , /a = 24 Nf eg gs- gs, RMniq, A) , (36) 

where Nf denotes a statistical fiavor factor which is equal to 1 or 2 for DHBs with two different or two identical heavy 
quarks. F^^, = {'i^qv — gtii/^) 75 is the Lorentz structure orthogonal to the photon momentum: q^ T-h^ — 0. Note that 
only the Lorentz structure F^^, survives in the HQL. As shown in Appendix A other possible structures vanish. The 
function Ri2{mg, A) reads 



d'^ki f d'^k2 2( \ niq + {k2 - ki)v 

''^' '''' ~ ' (2^ J ^2^ ^""^ {-kiv - ie)^{k2V - ie){ml - {k2 - fci)^ - leY 



RABimg,A)^ -— i- j:r-^<^'{z) , , ^. ,,,,,, ^.'^ :', '',, , ,, ,,,, (37) 



where z = — ^(fc^ — fcifc2 + A:|) and where A and B denote integer powers. The coupling constants gB* and (7^2 ^I'S 
given by 

-^ = -^ = 12Nf R2i{mg,A). (38) 

9bt ^Tb^ 



Using the Laplace transform 



c» 



$^(z)= / ds$^(s)e~"" (39) 



the integration over the virtual momenta fci and ^2 in i?As(wg, A) can be done. One obtains 

oo 

i?2i(m„A) = J^^^ J daida2da3— (^— + ^^ j $ (y) , 



oo 

^i^K.A) = (t|^/'^"i^"^^"3§(^ + ^^)$^(2/), (40) 



where 

^=4+"^' A^ = 3(,A^ ' i^ )- ^ ^ 

For the (| , 1) — > (^ ,1) transition the HQL hehcity amplitudes read 

^±iTi = -e*'^{Tl)u[v,±\)At.u-'[v,±^)^±M+M^ VW ^^ ' ^^^^ 



(43) 

The ratio of the HQL helicity amplitudes is given by ^±izri/^±i±i = ^VS, i-e. one has a pure Afl magnetic dipole 
transition. This coincides with the predictions of the nonrelativistic quark model (NQM) (see Appendix B), where 

if±wi = -V3i/±i±i =±^M+M_ ,/^. (44) 

±2=F1 ±2±1 ^3 ^ V Ml ^ ' 



The function /a appearing in P^ and (H51) is given by 



2^9 



/a-^/3, (45) 



/3^2m /,^-)l-^':; . (46) 



where 

Ri2{mq,A) 

'i?2l(TOq,A) 

In the HQL the helicity amplitudes, the form factors Fi and F2, and the decay rate of the (| , 1) -^ (^ ,1) transition 
read 



2u„ M2 

!=Fi ±2±i V3 \ Ml' 



H±^^i = -\/3iJ±i±i = ±-§/3M+Af_ 



^-^^-Wt' -— ^^fv/^. («) 



4 
3 
where 



Ti_,ls^;;K^ll^\ 



X = «M2M^4^. (48) 

a ~ 1/137 is the fine structure coupling constant and Mi, M2 are the masses of the parent and daughter baryon. The 
state i corresponds to baryon with symmetric heavy quark spin configuration. 

It is evident that our HQL rate results differ from the predictions of the NQM by the factor /3^ defined in ^^ . In 
the NQM one has (3=1 while in our covariant approach /? « 0.5. It is for this reason that our HQL predictions for 
the (| , 1) — > (5 , 1) + 7 decay widths are down by a factor of 4 compared to the predictions of the NQM. 
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IV. HYPERFINE MIXING AND RADIATIVE DECAYS OF MIXED STATES 

As mentioned in the introduction the origin of the hyperfine mixing in the double heavy baryons is the one-gluon 
exchange interaction between the hght and heavy quarks in the states containing two different heavy quarks — b and 
c. The one-gluon interaction leads to mixing of the states containing spin-0 and spin-1 heavy quark configurations. 
In this section we discuss in some detail the calculations of DHB radiative decays involving mixed states. We have 
three types of transitions: Bj^^ — > B^^, B^^ — >• B^^ and B^^ — > Bj^^. All three modes are quite interesting, because 
their study opens the opportunity to determine the mixing angle 6 and to measure the masses of the mixed states. In 
particular, the first mode B^^ — > B^^ is interesting since it is described by transitions between baryon components with 

the same spin configuration of the heavy quarks (i , 1) — )■ (i ,1) and (^ ,0) — )■ (^ ,0). Because now the photon 
can also couple to the light quark one will have a corresponding enhancement of the decay rates. The two other 
modes B^^ -^ B[^ and Bl^ -> B^^ involve mixing of the leading (|^, 1) -^ (i^, 1) and subleading (|^, 1) -^ (5^, 0) 
amplitudes and are therefore also important for an analysis of the mixing angle 0. The matrix elements for transitions 
involving mixed states are derived using the transition matrix elements of the unmixed states. In Appendix B we 
present the results of the NQM for transitions involving mixed states in terms of quark magnetic moments and the 
mixing angle 9. In our numerical calculations we differentiate between the mixing angles for the S-states and the 
fi-states using the predictions of the quark model pjj: 9^ ~ 25.5° and On ~ 25.9°. 

The last issue which we would like to discuss in this section is the HQL structure of transitions involving mixed states. 
As we have discussed in the previous section the leading contribution for the Bl^ — >■ B\^ and Bl^ — > B^^ transitions 
comes from the (| , 1) — ^ (i ,1) transition generated by the direct coupling of the photon with the hght quark (see 
diagram in Fig. 2(a)). The corresponding amplitudes are multiplied by the factor cos 9 for the Bl^ — >• B^^ mode and by 
— sin 6* for the B^^ — > B^^, mode. In the case of the B^^ — >■ B\^ transition the leading contribution is again generated 
by the direct light quark-photon coupling [Fig. 2(a)]. In this case one has to sum the two transitions involving a light 
quark spin-flip (i , 0) ^ (^ ,0) and (^ ,l)^(i ,1). The calculation of these leading matrix elements follows the 
treatment in the previous section. In particular, the leading contribution to the matrix element of the B^^ -^ B^^ 
transition is expressed though the same structure integral Ri2{'mq, A) as in the case of the (| , 1) — >■ (^ , 1) transition 

A^=^7^Ja, JA = 24 7V/e, sin(20B)^^±^i?i2(TO„A), (49) 

where gs and g'g are the coupling constants of the unmixed states (^ ,1) and (^ ,0): 

1 _ 1 
ffl' ^9b 
In the HQL the helicity amplitudes of the Bj^^ — > B^^ transition are given by 



2 -12 7V;i?2i(m„A). (50) 



i/±i±i =e*^(±l)fi(i;,±i)A>(z;,Ti) = Af+Af_ ^^ Ja . (51) 

After some straightforward algebra one can express the helicity amplitudes H±i±i in terms of the parameter /? derived 
in Eq. (gHI). One has 

^±i±i = ^ M« sin20 /? A/+ Af_ ./^ . (52) 

For /3 = 1 our helicity amplitudes coincide with the predictions of the NQM. We can also deduce the form factor F2 
[see the expression for the matrix element of the 1/2+ -^ 1/2+ transition (|A2[) ]: 

F2 = -'^ fi, sm29 P M, .^ . (53) 

Note that the form factor Fi defined in (jA2p vanishes due to gauge invariance. Finally the decay width for the 
^bc ~^ ^bc transition in the HQL reads 

T{B^,^Bl)^^Ktilsm^29, (54) 

which again coincides with the prediction of the NQM [see Eq. (jB12p ] when (3=1. 



involving mixed states i and ^ 
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Note that in the HQL our model also reproduces the model-independent results derived in m1^ for the decay rates 
t and i'. 

rj_h 11 ^ u^ sin 2Q , 

T3.1i ^ fj-l cos^e*, 

Fa ih ^ nl sin^^, 



r 



3_, 1 fi 



r^ 



tan^e . (55) 



V. RESULTS 

We now proceed to present our numerical results. We first present results on the radiative rates using finite heavy 
quark masses, i.e. we do not take the HQL for the matrix elements. Estimates for the decay widths are also given for 
the nonrelativistic quark model, in which, as described before, the wave functions have the same spin-flavor structure 
as our relativistic current considered in the nonrelativistic limit. Then we consider the HQL in both approaches. We 
choose the Gaussian form Eq. (|20|) for the correlation function of the double heavy baryons. Our results depend on 
the following set of parameters: the constituent quark masses and the size parameter A^. The parameters have been 
taken from a fit to the properties of light, single and double heavy baryons in previous analyses [15|: 

fnu{d) fns uric mb Kb ,rn\ 

0.42 0.57 1.7 5.2 2.5 - 3.5 GeV ' ^ ' 

All our analytical calculations have been done using the computer program FORM [27l |. 

In Table III we present detailed numerical results on the radiative rates of double heavy baryons using finite masses 
for the heavy quarks (exact results, second column) and in the HQL (third column). These are compared to the 
corresponding results of the NQM using finite masses for the heavy quarks (fourth column) . In column 5 we take the 
heavy quark limit of the NQM results by setting the heavy quark magnetic moments to zero. The dependence of our 
results on the size parameter Kb is indicated by error bars where the variation of A^ is given in Eq. (|56|) . Note that a 
smaller value of A gives bigger rates and vice versa. One can see that our finite heavy quark mass predictions are close 
to the results of the NQM. One has to keep in mind that the rate predictions are very sensitive to the mass difference 
AM — Ml ~ M2. In fact, one has F ^ (AM)'' (see Appendix B). The modes involving mixed states are enhanced 
by factors ~ 4 and ~ 2 in case of B^^ -^ B^^ and Bl^ — > B{^ transitions, respectively, while in case of Bl^ -^ B^^ 
transitions they are additionally suppressed by a factor ~ 10 due to reduction of the mass difference Mi — M2. In 
Table IV we present results for the rates involving mixed states in dependence on the mixing angle Ob varied from 
10° to 25". In Table V we compare our results with the results of a nonrelativistic quark model calculation [17|. The 
approach [17i] can be viewed as an extension of the naive NQM discussed before by taking into account baryon wave 
functions in configuration space. One should emphasize that the results of nonrelativistic quark models are in general 
frame-dependent. For example, the results can depend on whether one works in the parent baryon or daughter baryon 
rest frame. Also, it is difficult to maintain gauge invariance in nonrelativistic quark models. 

One final remark concerns the comparison of radiative and weak decays of DHBs. In [3] we have calculated the 
6 — >■ c semileptonic decays where we have shown that the corresponding decay widths are of the order of 10~^^ GeV. 
The radiative decay widths calculated in the present paper are much larger and lie in the range from lO^'* to 10^^ 
GeV. One would like to know how important the weak decays of DHBs induced by the qi -^ qj (d — > u or s — > u) 
light quark transitions are. For a precise analysis one would need to know the precise values of the masses of the 
DHBs including isospin-breaking corrections which are not available at present. Instead using the general formula for 
the semileptonic decay width one can obtain a rough estimate for the decay rates induced by light quark transitions 
where, for the sake of simplicity, we neglect the contribution of form factors, spin and flavor factors. For example, for 
the 1/2+ -^ 1/2+ transition one obtains (see e.g. [2S]) 

m^q,)^^m^AM\ (57) 

where Gp = 1.16634 x 10~^ GeV~^ is the Fermi constant, Vckm is Gabibbo-Kobayashi-Maskawa (CKM) matrix 
element (|KidP — 1 and |T4sP = 0.051), and AM = Mi — M2 is the difference of the masses of initial and final 
baryons. We know from data on the mass differences of light and heavy-light baryons that, approximately, the mass 
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difference AM does not exceed the mass difference of the corresponding hght quarks in these baryons. Therefore, in 
the expression for T{qi — >■ Qj) we substitute AM < m^ — rriu for d ^ u transitions and AM < rUs — m^ for s ^ u 
transitions and considering this to be an approximation of the upper hmit for the corresponding decay rates. Using 
upper hmits m^ — 77i„ < 10 MeV and nis — rriu < 200 MeV we obtain: 

r(d -^u) < 10^22 QgY ^ Y{s -^ u) < 10"^^ GeV . (58) 

Based on this rough estimate one concludes that the semileptonic decays of DHBs induced by the hght quark transitions 
d — >■ u and s — >■ u are suppressed by more than 4 orders of magnitude in comparison to their 6 — >■ c counterparts and 
are even more suppressed (more than 10 orders) compared to their radiative decays. 

VI. SUMMARY 

We have analyzed the radiative decays of double heavy baryons using a manifestly Lorcntz covariant and gauge 
invariant constituent quark model approach. Our main results can be summarized as follows. We have derived results 
for the radiative transition matrix elements of double heavy baryons for finite values of the heavy quark/baryon masses 
and also in the HQL limit of infinitely heavy quark masses. We have discussed in detail radiative transitions involving 
DHB states subject to hyperfine mixing. We have presented an extensive numerical analysis of the decay rates for 
finite masses and in the HQL limit including numerical results on mixing effects. Our results were compared with 
the predictions of a nonrelativistic quark model including again hyperfine mixing effects. We find that the inclusion 
of hyperfine mixing effects has a profound influence on the pattern of radiative decays of DHBs. Since the calculated 
rates depend very sensitively on the exact mass values of the mixed and unmixed DHB states (F ~ (Mi — M2)^) one 
must wait for an accurate determination of the masses of the DHB states before one can extract information on the 
mixing angles from the decay data. 
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Appendix A: Spin— kinematics of radiative decays 

In this Appendix we write down covariant expansions for the current-induced electromagnetic transitions involving 
the (1/2+) and (3/2+) baryon states. We thereby define sets of invariant vector transition form factors. We then 
define helicity amplitudes which are expressed in terms of linear combinations of the invariant form factors. One of 
the advantages of using helicity amplitudes is that one obtains very compact expressions for the decay rates (see e.g. 
[28|,|29|). In addition, the helicity amplitudes contain the complete spin information of the process and are thus well 
suited for the computation of spin observables. 

In the radiative decays of double-heavy baryons the momenta and masses are denoted by 

Bi(pi,Mi)^S2(p2,M2)+7('z) (Al) 

where pi — P2 + q- For the invariant form factor expansion of the 1/2+ — > 1/2+ matrix elements of the vector current 
J^ one obtains 

Transition i — > i : 

M^ = (i?2|J^|Si) - u{p2,S2)U^F^{q^)-ia^,^F2{q^) + ^F^(q^)\u{p^,si) . (A2) 

L ivl\ 1V1\ J 

Similarly one has 
Transition § — J' 5 : 

M, = {B2\J,\Bl) = u{p2^s2)[9o.,F,{q^)+l/-^j2{q^) + ^^F,{q') + ^J^4(g')]75Zi"(pi, ^i) , (A3) 
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where a^u — {i/2){j^"ii, — Jujfj.) and all 7 matrices are defined as in Bjorken-Drell. One should emphasize that the 
above invariant form factors are constrained by gauge invariance relations (see e.g. the detailed discussion in [30]). 

Next we express the vector helicity amplitudes H\^\ in terms of the invariant form factors Fi, where A^ = ±1 
and A2 — ±1/2, ±3/2 are the helicity components of the on-shell photon and the daughter baryon, respectively. The 
pertinent relation is 

Hx,x^^M^{X2)e*^'{X^). (A4) 

Angular momentum conservation fixes the helicity Ai of the parent baryon according to Ai = A2 — A^. We shall 
work in the rest frame of the parent baryon Bi with the daughter baryon B2 moving in the positive z-direction 
such that p^^ = iMi,0), pt^ = (S2,0,0,|p2|) and q^' = (go, 0,0, -|p2|), where qo = IP2I = (Mf ~ M|)/(2Afi) and 
E2 = Ml - go = (Mf + M|)/(2Mi). 
The J = ■^ baryon spinors are given by 

ui{pi,±l) = v/2Mr(j^) (A5) 

where x+ ^ ( n ) ^^^^ 't- = ( 1 ) ^'"^ two-component Pauli spinors. 
The J — ^ baryon spinors are defined by 

u^{p,s*) ^^{lX-s\-s*)ef,{p, X)u{p, s) . (A6) 

They satisfy the Rarita-Schwinger conditions 

j^u,,{p,s*)^p''u^ip,s*)^0, (A7) 

where (lA^sJls*) is the requisite Clebsch-Gordan coefficient, e^{p,X) is the spin 1 polarization vector and u{p,s) are 
the usual J — ^ spinors defined above. In particular, the J = | spinors with helicities A — ±3/2, ±1/2 read: 

U,.{P,±^) = ^J-ef,ip,0)u{p,±-j +^-e^{p,±l)u(p,T^) 

The polarization vectors corresponding to the parent and daughter J — ^ baryons are given by: 

e^(pi,0) = (0,0,0,1), e^(Pi,±l) = -i=(0,Tl, -*, 0) , 

e*^(P2,0) = ^(|p2|,0,0,ii;2), e*'^(p2,±l) = -^(0,Tl,«,0). (A9) 

The polarization vectors of the on-shell photon read 

e-*'^(±l) = -^(0,±l,z,0). (AlO) 

where the "bar" on the polarization vector denotes the fact that the photon is moving in the negative z-direction. 
The polarization vectors satisfy the Lorentz condition 

g^6-*''(±l)=0. (All) 

Using above formulas for the spin wave functions with definite helicities one can then express the helicity amplitudes 
Hx^x through the invariant form factors by calculating Hx^x — M^{X2)e*^{X^). One obtains 

H,.,,^-F2V2^, (A12) 



(A8) 
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where M± = Mi ± M2. 
Transition f ^" ^ 



1 / Ml 



V3 V ' M1M2' 

i/±i^i = ±M_^i. (A13) 

Note that one has an exphcit factor of M_ in all the helicity amplitudes which corresponds to the appropriate p-wave 
threshold factor |p| ^ M-. The decay width is given by 

_ a Mf-Mj 
where the 'Hsi^s2 ^-rc bilinear combinations of the helicity amplitudes: 

2 2 2 2 

^3^1 = \Hu\^ + \H_.A^ + \H^A'' + \H_u\\ (A15) 



The overall dependence of the rate on the mass difference M_ can be seen to be given by |p| ^'+^ ^ M^ for 1 = 1. In 
the HQL the baryon spinors simplify. For example, the HQL J = 5 baryon spinors read 

(«,±i)=y2M2(xL,0), u^(v,±^=^2mAiA . (A16) 



U2 

For the | ~^ 5 transition in the HQL one has four possible Dirac strings in the matrix elements which are 75 , 75 ^, 
75 t*^^!) ^nd 75 4 t*{^'<)- It is easy to show that the two Dirac strings 75 and 75 ^ ^*{\j) vanish when sandwiched 
between the HQL baryon spinors: 



; 2) (75,75^ ^*(A^) ju^(^w,±-j =U2{y,±-j[-i5,-i^i^*{\^)ju^,[y,^- 



The remaining strings 75 ^ and 75 ^* one has 



U2(«,±^)75 iu^(y,±^ =tM+M- ^^e^(^;,±i), 

W2(w,±-)75 fiu^,(^v,T^) =T\j^M+M^ y'_ie^(„,=pl), 

"2(w,±-)75 ^*(t1)u^(w,±-) =0, 

U2(«,±J)75 r(±l)"M(«-T^) =±4^ V^iM2e^(«,0). (A18) 



In the calculation of the helicity amplitudes iJ_j_i_|_]^ and i^±i^i one can make use of the HQL identities 



^2(«,±i)75 r{±l)qW{v.T\) ^±^M+M^ ^^ (A19) 

Appendix B: Nonrelativistic quark model: spin— flavor wave functions, radiative decay constants and widths 

of double heavy baryons 

In this Appendix we present results on the radiative decay amplitudes and widths of the DHBs in the nonrelativistic 
quark model. As emphasized before the nonrelativistic quark model is based on the spin-flavor wave functions which 
arise in the nonrelativistic limit of the relativistically covariant double-heavy three-quark currents with quantum 
numbers J^ = -k and | . The corresponding quark model spin-flavor wave functions are given in Table 2, where 
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we use the following notation for the antisymmetric XaW ^^'^ symmetric XsW^ X*!-^) spin wave functions where A 
is the helicity of the baryon state: 

X.© - /|{t(n-;t)}, x.(4) - /f {;(n-;t)}, 

^■'Q ^ V^{2in-t(n + ;t)}, x.(-^) - -y||2m-;(n + ;t)}, (bi) 

x;{l) = ]llUn + m + m\,, ^^(-^) = y^ In; + ;n + utj . 

In (JB1|1 we use the ordering {qQiQ2}- Next we relate the DHB radiative decay amplitudes to the nonrelativistic 
amplitudes G^^Xj, where A^ = ±1 and A2 = ±1/2, ±3/2 are the helicity components of the on-shell photon and the 
daughter baryon, respectively (see details in Appendix A). In order to evaluate the spin flip matrix elements between 
the baryon states we make use of the spin-flip (spin raising/lowering) operator 

3 

5±p = -y2^A*,(a±),, (B2) 

1=1 

where i runs over the three constituent quarks. The spin flip matrix elements are given by 

Ga.ti = (52(A2)|5±p|Bi(Ai)), (B3) 

where Ai = A2 — A^ is the helicity of the parent baryon and cr± = cti ± ia2- fJ-i — Ci/ [2mi) is the i-th quark magnetic 
moment and e^ and rrii are its charge and mass, respectively. 
For the amplitude Gx^x one obtains 



Transition i " ^ 



2 ~^ 2 



Transition | 



Transition 



2 ^^ 2 



S 



Transition ^ ~^ 2 



Transition | -^ 5 



Transition 



3^1 

22 



h 



G±^-±i = -\l^{f^c-f^b)- (B4) 



G'±i^i = -V3G±i±i = ±(Aic-M6). (B5) 



G±.^r = -^/3G±.±, = ±-| (m, - ^^^1±^) . (B6) 



9 /o / /^ 

G±i±i = -IT- sm29B ( Mg - ^ (a*c - ^J■b) cot29B ] ■ (B7) 



G±izpi = -\/3G±i±i =±^ cose's Uq- Ytan+-y tan_ ) . (B8) 



G±izpi = -\/3G±i±i =±-^ sine's Ug-yCot_-yCot+) , (B9) 
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where tan± = 1 ± tan^^-v/S and cot± = 1 ± cot ObVS- 

The states 5 : 2 ' 2 ^^'^ 2 correspond to baryons with antisymmetric, symmetric and mixed heavy quark 
spin configuration, respectively. In particular, i and ^ correspond to the mixed states Bfi and _B;. The index q 
corresponds to the light u,d or s quark. 

In terms of the nonrelativistic amplitude Gx^x the helicity amplitudes H\^\ (sec details in Appendix A) for the 

^ ^ radiative transitions are defined by: 



i and 2 



i?A,A, -|P2|V^^G. 



A2 A^ 



(BIO) 



where Ni — \/2Mi is the extra factor acquired in the nonrelativistic quark model due to different normalization of 
states in the relativistic and the nonrelativistic theory. Therefore, we have 



M2 
V -'^i 

The radiative decay widths for the four possible si ^- S2 spin transitions are given by: 

2 



FlA iS 



r ih u 



2^2 v ^^ 



-Kill sm^2eB (l-^ ^^^-^ cot 29b 



r3^is = -if/ia 1 



2flq 



(Bll) 



where 



Ts^ii = -K ni cos^Ob 1 - -^-^tan+--^— tan_ 
3 ^ V 2/Xg 2^iq 



3_i 1 '' 
2 ~^ 2 



\K^,ls\n^eB(l-^cot^-^coU 

3 ^ \ 2fJ,q 2^5 



2\3 



K = aM2 



{Ml - M l) 
6Mf 



and a ~ 1/137 is the fine-structure constant. 

It is evident that the widths of the subleading processes i' 
{mq/mcf compared to the widths of the leading process | -^ -. 



(B12) 



(B13) 



■i are suppressed by a factor of 



Appendix C: Mass operator and radiative vertex functions of double heavy baryons 

The baryon mass operators S^ (p) and E^'C (p) are given by 

Sb(p) = 67V/ fdki23<^\zo)Rs{k+,k+,k+), 
i?s(ri,r2,r3) = ri/S'g(r3)rii tr T2fSQ^{r2)T2iSQ^{-ri) 



(CI) 



S^-: (p) ^ 6 I dk,23 <^'{zo) R'^'{kt,k+, k+) , 
R'^ (ri,r2,r:i) = ri/S'5(r3)rH tr T'^fSQ^{r2)V2,SQ,(-ri) 



(C2) 
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Here and in the following Nf denotes a statistical flavor factor, which is equal to 1 or 2 for DHBs with two different 
or two identical heavy quarks, respectively, and T = 7''r'l'7''. We have introduced the abbreviations 



and 



dki23 = 



d^kid^k2d'^k: 



3 .4 



{2n)»i 



8,-2 



r(fci + fc2 + fcs) , zo = -- {kf + k^ + ki) 



kf = ki+pwi, fc^+ = ki+p'wi, Li = -[kj -/^k_ 



iW, 



j=i 



zi{q) = Liq- -q^(wl+W2W3+wl). 



Mq) = L2q--q (W1+W1W3+W3) 

2 
zsiq) = L^q- -q'^{wl+WiW2+wl) 



Ri^,Ai{ri,r2,r3,q) = -TifSq{r3)Tiitr T2fSQ^{r2)T2tSQj-ri)'j^SQ^{-ri+q) 



-Rm,A2(^1'^2,?-3,9) = TifSg{r3)rutr r2fSQ^{r2 ~ qh^SQ^{r2)T2iSQ,{-ri) 
Ra,A3{ri,r2,r3,q) = TifSqir^ - q)-f^Sq{r3)rutr T2fSQ^{r2)T2iSQ^{-ri) 



^^,.,Ai(^i>^2,?-3,<z) = -ri/5'g(r3)rHtr T2fSQ^{r2)T2i.,.SQ,{-ri)j-^SQ,{-ri + q) 



(C3) 



^Mi.,A2(''i'''2,''3,'Z) == ri/^q(r3)rHtr r2fSQ^{r2 - qh^SQ^{r2)T2i.iySQ^{-ri) 



^A^,.,A3(^i'^2,r-3,g) = TifSqirs ~ q)'^^Sq{r3)rutr T2fSQ^{r2)T2i^t.SQ,{-ri) 



Ru,T.{ri,r2,r3) = ri/^q(r3)rH tr T2fSQ^{r2)T2i,ySQ^{-ri) 



(C4) 



1 + 1 + 
In the following we present explicit expressions for the electromagnetic vertex function. In case of the (2 , 0) — > (^ ,1) 

and (| ,1) — > (2 ,0) transitions the expressions for the nonvanishing contribution of the triangle diagrams in 

Figs.2(a)-2(c) (terms A^ ^ and A^^ ^, respectively) read 

^iA{PuP2) = 6NfgB,9B2 dki23 ^ Ci $(2:0) $(zo + z,(g)) i?^^A^(fc+, fc+, fc+, g) , (C5a) 

f ^ 

^tf,A{Pi^P2) ^ QNfgBigB2 dki23 ^ e^<^{zo)<^(zo + z,{q)j Rj^^^^^{kf,kf,k^,q), (C5b) 



flu 



where ei = cqj, 62 = eg^ and 63 — eg. For the (| , 1) — > (i , 1) transition the electromagnetic vertex function A 
obtains contributions from the triangle diagram A^^. a [Fig.2(a)-2(c)], the left and right bubble diagrams A^^^^^^^^ 
[Fig. 2(d)] and A^^^ ,^^,^^ [Fig. 2(e)], and the pole diagram Aj^^ j [Fig. 2(f)]. The corresponding contributions are given 



by 
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Kl^.huhr^{Pl:P2) 



-%NfgBi9B, f dki23 ^e»i4$(zo) j dt^' (z^ + tz,[-q)^ R^^^{k' + ,k' + ,k'^+) , {C&&) 

3 1 
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FIG. 1: Diagram describing the double heavy baryon mass operator. 




FIG. 2: Diagrams contributing to the radiative transitions of double heavy baryons. 
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Table I. Classification and mass values of double-heavy baryons. Mass values are based on [13 
except for the Hcc mass which is taken from 23]. 



Notation 


Content 


r 


Sd 


Mass (GeV) 


^cc 


q{cc} 


1/2+ 




3.5189 


Si,c 


q{bc} 


1/2+ 




6.933 


^bc 


q[hc] 


1/2+ 





6.963 


'^bb 


q{hb} 


1/2+ 




10.202 




q{cc} 


3/2+ 




3.727 


^bc 


q{bc} 


3/2+ 




6.980 


—66 


q{hb} 


3/2+ 




10.237 


n,c 


a{cc} 


1/2+ 




3.778 


^bc 


s{bc} 


1/2+ 




7.088 


fiL 


s[bc] 


1/2+ 





7.116 


^66 


s{bh} 


1/2+ 




10.359 


f^cc 


s{cc} 


3/2+ 




3.872 


^*bc 


8{bc} 


3/2+ 




7.130 


^Ib 


s{bb} 


3/2+ 




10.389 



Table II. DHB wave functions. 



Baryon 


Wave function 


Baryon 


Wave function 


^cc 


ICC XsW 


Occ 


sec XsW 


'S.bb 


qbb Xs(A) 


Qbb 


sbb XsW 


Sfcc 


^^q{bc + cb) Xs(A) 


nt. 


^s{bc + cb) XsW) 


—be 


^q{bc-cb) x^(A) 


^L 


^^s{bc-cb) xA^) 


^cc 


-qcc x*sW 


ficc 


-sec x*sW 


^bb 


-qbb x*sW 


^Ib 


-sbb xlW 


^bc 


-^q{bc + cb) xlW 


^L 


-^s{bc + cb) x*sW 



Table III. Radiative decay widths of DHBs in keV. 
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Decay mode 


Exact results 


HQL 


NQM 


NQM + HQL 


"6c ^ "6c 


(1.56 ± 0.08) X 10"^ 





1.35 X 10~2 





"6c ^^ "6c 


(1.56 ± 0.08) X lO"'' 





1.35 X 10~2 





fiL -^ ^bc 


(1.26 ± 0.05) X 10"^ 





1.10 X 10^2 





-be "^ "6c 


0.14 ± 0.03 


~ 0.17 


0.34 


0.53 


who _, ^!0 
—be ^^ "6c 


0.31 ± 0.04 


~ 0.04 


0.26 


0.13 


^t -^ «L 


0.21 ± 0.02 


-0.02 


0.15 


0.06 


"6c ^^ "6c 


(0.28 ± 0.01) X lO"'' 





0.25 X 10-2 





"6c -^ —be 


(0.28 ± 0.01) X 10^2 





0.25 X 10^2 





^bc -> «L 


(0.16 ± 0.01) X 10"^ 





0.14 X 10^2 







23.46 ± 3.33 


20.53 ± 0.79 


36.22 


63.88 


^-"ce ^ '— 'cc 


28.79 ± 2.51 


5.13 ± 0.20 


35.65 


15.97 


n*, ^ rice 


2.11 ± 0.11 


-0.29 


2.42 


0.87 


"6c ~^ -be 


0.49 ± 0.09 


~ 0.27 


0.67 


0.83 


"6c —^ "6c 


0.24 ± 0.04 


~ 0.07 


0.30 


0.21 


fibc ^ fibc 


0.12 ± 0.02 


~ 0.03 


0.13 


0.08 


"6c ^ "6c 


0.46 ± 0.10 


-0.37 


0.69 


1.14 


"6c ^ "6c 


(0.15 ± 0.02) X 10-2 


~ 0.03 X 10^2 


0.16 X 10^2 


0.08 X 10"2 


"6c ^ "6c 


0.51 ± 0.06 


~ 0.10 


0.59 


0.28 


"6c ^^ "6c 


(0.02 ± 0.02) X lO""* 


~ 0.06 X 10"^ 


0.01 X 10"* 


0.19 X 10"* 


ni ^ nL 


(0.29 ± 0.03) 


~ 0.03 


0.30 


0.12 


fisc ^ ^t 


(0.01 ± 0.01) X 10-* 


~ 0.01 X 10"^ 


0.01 X 10"'' 


0.03 X 10"* 


—bb ^^ "66 


0.31 ± 0.06 


~ 0.11 


0.38 


0.35 


"66 ^^ "66 


(5.87 ± 1.42) X 10^2 


~ 2.8 X 10^2 


7.34 X 10"2 


8.69 X 10"2 


^66 -^ ^66 


(2.26 ± 0.45) X 10^2 


~ 1.0 X 10^2 


2.36 X 10"2 


2.97 X 10"2 
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Table IV. 9b dependence of radiative decay widths involving mixed DHBs in eV. 



9b 


Decay mode 


Exact results 


HQL 


NQM 


NQM + HQL 


10" 




0.2 ± 0.2 


- 10 


3 


31 




170 ± 15 


-3 


43 


8 


^t ^ ^L 


130 ± 10 


- 1 


28 


3 


"be ^ "6c 


0.07 ± 0.01 


- 0.3 


5 


1 


^»0 , "hO 
"be — ^ "6c 


-0.1 


- 0.1 


0.7 


0.3 


fiftc ^ ^t 


-0.03 


-0.02 


0.5 


0.1 


"be ^^ "be 


687 ± 132 


276 ± 11 


626 


859 


^»0 , —H) 
"be ^^ "be 


475 ± 65 


69 ±3 


360 


215 


f7L ^ nl 


263 ± 26 


- 28 


164 


84 


15° 


"be ~^ "6c 


16 ± 7 


27 ± 1 


27 


84 


"6c ~^ —be 


224 ± 24 


— 7 


73 


21 


fifcc ^ f^6c 


166 ± 14 


-3 


44 


9 


"be ^^ "6c 


0.9 ± 0.1 


- 0.6 


6 


2 


^»0 , -^hO 
"be ^^ "6c 


-0.05 


-0.1 


0.2 


0.5 


f^6c ^ ^L 


-0.02 


-0.1 


0.2 


0.1 


"be ^ "be 


619 ± 122 


289 ± 12 


619 


900 


"be ~^ "be 


493 ± 66 


72 ± 3 


404 


225 


f^L ^ "be 


276 ± 26 


- 29 


187 


88 


20" 


"be ^ "6c 


63 ± 19 


64 ± 2 


100 


200 


"6c -^ "6c 


272 ± 31 


16 ± 1 


126 


51 


nL -^ nL 


195 ± 17 


- 7.3 


72 


22 


"be ~^ —be 


1.2 ± 0.1 


- 1.5 


5 


2 


—be — ^ "6c 


- 0.01 


- 0.2 


0.01 


0.5 


fifcc ^ ^be 


- 0.01 


-0.03 


0.03 


0.1 


"be ~^ —be 


546 ± 111 


313 ± 12 


632 


973 


-be — ^ —be 


504 ± 65 


78 ± 3 


468 


243 


^te ^ ^L 


286 ± 27 


- 32 


223 


96 


25" 


= '1+ ^ ■=' + 
"be ^ "6c 


135 ± 34 


91 ± 4 


177 


284 


"6c -^ —be 


306 ± 37 


223 ± 1 


142 


71 


^t -> nL 


213 ± 20 


- 18 


131 


52 


"be ^^ "6c 


1.4 ± 0.2 


-0.3 


2 


1 


^*0 _, ■^hO 
^be ' ""be 


0.001 ± 0.001 


- 0.08 


0.004 


0.2 


f^6c ^ «6c 


0.001 ± 0.001 


- 0.02 


0.002 


0.04 


"be ^^ "be 


469 ± 99 


373 ± 15 


680 


1115 


^»0 , "«) 
"be ^^ "be 


507 ± 64 


90 ± 3 


578 


278 


f7L ^ nl 


291 ± 26 


- 37 


281 


HI 
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Table V. Radiative decay widths of DHBs in keV. Comparison with the quark model [17 



Decay mode 


Quark model [17] 


Our results 


"6c ^ "6c 


0.992 X 10"^ 


(1.56 ± 0.08) X 10-2 


"6c -^ ~^bc 


0.992 X 10-2 


(1.56 ± 0.08) X 10-2 


fifcc ^ ^b. 


3.69 X 10^2 


(1.26 ± 0.05) X 10-2 


"6c "*■ "6c 


12.4 X 10-2 


(14 ±4) X 10-2 


^bc -^ ^bc 


20.9 X 10-2 


(31 ± 4) X 10-2 


^t ^ "L 


8.52 X 10-2 


(21± 2) X 10-2 


"6c -^ "6c 


4.04 X 10-2 


(0.28 ± 0.01) X 10-2 


--be -^ ^bc 


4.04 X 10-2 


(0.28 ± 0.01) X 10-2 


fiftc ^ "L 


3.69 X 10-2 


(0.16 ± 0.01) X 10-2 


"6c ^ "6c 


1.05 


0.49 ± 0.09 


^»0 , -^0 
"6c ^^ "6c 


0.505 


0.24 ± 0.04 


0,1c ~> ^6c 


0.209 


0.12 ± 0.02 


"6c ^^ "6c 


0.739 


0.46 ± 0.10 


--be ~^ "be 


6.05 X 10-2 


(0.15 ± 0.02) X 10-2 


"6c ^^ "6c 


1.03 


0.51 ± 0.06 


^»0 , —hO 
"be ~* "6c 


0.12 X 10-2 


(0.02 ± 0.02) X 10-"* 


f^6c -^ ^ie 


0.502 


0.29 ± 0.03 


ni ^ nt 


0.31 X 10-2 


(0.01 ± 0.01) X 10-* 



